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ABSTRACT 

The first complete orbital solution for the double-lined spectroscopic binary system 
7^ Velorum, obtained from measurements with the Sydney University Stellar Interfer- 
ometer (SUSI), is presented. This system contains the closest example of a Wolf-Rayet 
star and the promise of full characterisation of the basic properties of this exotic high- 
mass system has subjected it to intense study as an archetype for its class. In combi- 
nation with the latest radial-velocity results, our orbital solution prod uces a distance 
of 336I7 pc, significantly inore d istant than the Hipparcos estimation (jSchaerer et al.l 
[Hill van der Hucht et allllOOTl ). The ability to fully specify the orbital parameters 
has enabled us to significantl y reduce uncertaint ies and our result is consistent with 
the VLTI observational point (iMillour et al.l l2006'). but not with their derived distance. 
Our new distance, which is an order of magnitude more precise than prior work, de- 
mands critical reassessment of all distance-dependent fundamental parameters of this 
important system. In particular, membership of the Vela 0B2 association has been 
reestablished, and the ag e and distan c e are also in good accord with the population of 
young stars reported by iPozzo et a l. ( 2000) . We determine the 0-star primary com- 
ponent parameters to be My(0) = -5.63±0.10mag, R(0) = 17±2R0 and A^(0) = 
28.5 ± 1.1 Mq. These values are consistent with calibrations found in the literature if 
a luminosity class of II-III is adopted. The parameters of the Wolf-Rayet component 
are M„(WR) = -4.33 ± O.lTmag and 7W(WR) = 9.0 ± 0.6Mq. 

Key vi^ords: stars: individual: 7^ Vel - stars: fundamental parameters - stars: Wolf- 
Rayet - binaries: spectroscopic - techniques: interferometric 



1 INTRODUCTION 

The high-luminosity Wolf-Rayet (WR) stars are charac- 
terised by bright emi ssion line spectra pro duced in hot, high- 
speed stellar winds Ijvan der Huchtl ll99^). Thought to em- 
body the final stable stage in the evolution of very massive 
stars, t hey are candidate p recursors to Type Ib/Ic super- 
novae (Ivan der Huchtl bOOll ). The fast {voo lO^kms"^) 
stellar winds of WR stars are responsible for prodigious 
mass loss (10~* Mq yr~^), eventually stripping the star of 
its outer layers and forming their characteristic He-rich 
emission-line spectrum. Despite the ephemeral nature of this 
phase, these stars are important not only for understanding 
the evolution of massive stars but also as tracers o f galactic 
structure and star formation (Ivan der Huchtll200ll ). 

The double-lined spectroscopic binary 7^ Velorum 
(HR 3207, HD 68273, WR llQ) contains the brightest 
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example of the WR type in the sky and has a com- 
bined visual magnitude V ~ 1.8 with an 0-star compan- 
ion. It has a Hipparcos deter r nined d istance of 258t.t] PC 
dSchaerer. Schmutz fc GrenonI 1 19971 : Ivan der Hucht et al.l 
Il997f) . placing 7^ Vel closer to Earth t han any other WR 
star by a factor of a pproximately two (|Morris et al.ll2000l : 
Setia Gunawan et al.l 120011 ) . Its proximity is not the only 
reason that it is a key target for astronomers: the orbital 
motion of the two stars makes 7^ Vel one of the few systems 
for which a direct mass determination of the components 
can be achieved. As such, 7^ Vel provides critical tests for 
the theory of WR and massive star evolution. It is therefore 
not surprising that numerous studies of this system have 
been made, but in spite of its seemingly textbook status, 
considerable uncertainties remain. 

Spectroscopic analysis of 7^ Vel has been hampered 
by entanglement of the 0-star and WR wind emission: 
all absorption features of the brighter O-star compo- 
nent are blended with emission lines from the WR wind 
(|De Marco et all I2OO0I V Accounting for the WR emission 
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is not trivial because a) multiple emission lines may con- 
tribute to the blend, b) the shape of the emission line 
may differ significantly from a Gaussian, c) the emis- 
sion lines vary in time with some excess emission at- 
tributed to wind-wind collision phenomena, d) the ex- 
tended wi ngs of the strong O-star absorption are lost in 
the blend llSchmutz et al.ll 19971: iDe Marco fc Schmut j | l999l : 



Table 1. Adopted parameters of reference stars used during ob- 
servations. 



lEversberg. Moffat fc Marchenkolll999l '). This has resulted in 
numerous luminosity and spectra l class identificati ons of 
the O-star in the literature: Q7. 5 llGanesh fc Bappr] 
08 (iBaschek fc Scholzlll97ll) . 09 I JConti fc Smithf 



1967 ). 



1972 V, 



Q8II I l|Schaerer et all! 19971 ) and 07.5 III-V (|van der Huchi] 
I2OOII ) .The WR compo nent ha s retain ed its original classifi- 
cation of type WC8 bv lSmithI |l968i). 

The first to determine the period and some of the 
orbital elements by ra dial velocity measurements were 
iGanesh fc Bappul l| 19671). Further spectroscopic studies, for 
examp le lNiemela fc Saha.del (ll980D . |Pike. Stickland fc Willij 



and Moffat et al. ~(|l986l ). continued to add to 



1198 

knowledge of the system but discrepancies could not be re- 
solved and observations from the lUE co uld 'only be seen as 
mudd ying some already murky waters' (jStickland fc Llovdl 
I1990I ). By correct ing the O-star absorption lines due to 
the WR emission, ISchmutz et al.l l| 19971 ) believe they have 
the most definitive spectroscopic orbit to date. However, 
this type of analysis cannot determine all of the orbital 
elements. For example, the incli nation has been estim ated 
from the polarimetric analysis of lSt.-Louis et all lll987l) an d 



further constraints we r e imp osed b v ISchmutz et al, 
iDe Marco fc Schmutzl (|l999h and iDe Marco et al 



who used the Hipparcos distance, stellar models and spec- 
troscopy. 

Due to the high angular resolution needed, only in- 
terferometric techniques can currently be employed to ac- 
curately constrain the remaining orbital parameters. How- 
ever, the southerly declination (— 47°20') of ■y'^ Vel places 
it out of re ach for many interferometers of sufficient resolv- 
ing power. iHanburv Brown et al.l (|l970l ) used the Narrabri 
Stellar Intensity Interferometer (NSII) to study this sys- 
tem, determining the angular semi-major axis to be 4.3 ± 
0.5 mas. Unfortunately, the signal-to-noise of their obser- 
vations necessitated th e adoption of the parameters of 
iGanesh fc Bappul (Il967l) which have now been superseded 
by those of Schmutz et al.l ([l997'). Furthermore, due to the 
lack of observations, the inclination could not be adequately 
constrained and a value i — 70° was assumed from theo- 
retical arguments. The Very Large Telescope Interferometer 
(VLTI) has observed 7^ Vel a nd produced a single angular 
separation and position angle (|Millour et al.ll200q ) but has 
yet to improve the orbital parameters. 

Reliance on distance estimations and stellar models has 
paved the way for further analysis of fundamental parame- 
ters, although this path has been fraught with controversy. 
The pie-Hipparcos distance of 7^ Vel was inferred from the 
assumed membership of the Vela OB2 association or from 
statistical absolut e mag nitudes of similar or nearby stars. 
iBaschek fc Schol3 (Il97ll ) ass umed equ idistance with the vi- 
sual companion 7^ Vel while IConti fc S mith ( 1_97^) adopted 
460 pc from estimates of nearby s t ars. T he interferomet- 
ric study of iHanburv Brown et al.l (Il970l) c ombined with 
the spectroscopy of Ganesh fc Bappul 1 19671 ) produced the 
first dynamical parallax of this system yielding 350 ± 50 pc, 



HR 


Name 


Spectral 


V 


UD Diameter 


Separation 






Type 




(mas) 


from 7^ Vel 


3090 


J Pup 


B0.5 lb 


4.24 


0.18 ±0.03" 


2.96° 


3117 


X Car 


B3I Vp 


3.47 


0.36 ±0.06" 


6.02° 


3165 


C Pup 


05f 


2.25 


0.41 ±0.03'' 


7.41° 


3468 


a Pyx 


B1.5 III 


3.68 


0.27 ±0.06" 


15.54° 



Notes: " Based on NSII measurements of similar type stars. 
Hanbury Brown. Davis fc Allen (1974) . 



although this distance is affected by the assumed values 
of some orbital parameters. In the 1980s, the distance to 
7^^ Vel was assumed to be equal to that of the Vela OB2 
association (450pc) with sufficient certainty as to be used 
in the determination of ga lactic WR intrinsic parameters 
(|van der Hucht et al.|[l98g ). The publication of the signif- 
ica ntly different Hipparc os pa rallax distance of 25 8 l^j pc 
by ISchaerer et~aLl l|l997h and Ivan der Hucht et all (| 19971 ) 
caused a critical reassessment of all 7^^ Vel parameters that 
are distance dependent - the absolute magnitude of the 
components, the spectral class of the O-star primary, the 
mass-loss of the WR star etc. The Hipparcos distance has 
been challenged by the discovery of a low-mass, pre-main- 
sequence stellar associatio n in the dire c tion o f 7^^ Vel and 
the Vela 0B2 association. IPozzo et al.l (I2OO0I ) argue these 
low-mass stars are associated with 7^ Vel and approximately 
coeval at a distance consistent with earlier estimates of 360- 
490 pc. The single measurement by the VLTI has yielded a 
measure of the distance to 7^^ Vel of 3681^3 pc (|Malbet et al.l 
2006; Millour et al. 2006) placing the system farther from 
the Hipparcos distance and closer to the Vela OB2 associa- 
tion. 

In Section [3] we provide the first complete orbital solu- 
tion for 7^^ Vel determined by long-baseline optical interfer- 
ometry. In combination with the latest radial velocity mea- 
surements, the tightly constrained distance, spatial scale and 
mass of the components are quantified and compared to pre- 
vious estimates in Section |4] The details of our observations 
and a description of the parameter fitting procedure are de- 
scribed in Section [2] and Section Irrespectively. 



2 OBSERVATIONS AND DATA REDUCTION 

Measurements of the squared visibility (i.e. the squared 
modulus of the normalised complex visibility) or were 
completed on a total of 24 nights using the Sydney Univer- 
sity Stellar Interferometer (SUSI. lOavis et al]|l999l ). Inter- 
ference fringes were recorded with the red beam-combining 
system using a filter with centre wavelength and full-width 
half-maximum of 70 nm and 80 nm resp ectively. This sys- 
tem was outlined in iTuthill et all (|2004l ) and is to be de- 
scribed in greater detail in Davis et al. (in preparation). 

Interference fringes produced by a pupil-plane beam- 
combiner were modulated by repeatedly scanning the optical 
delay about the white light fringe position. The two outputs 
of the beam-combiner were detected by avalanche photo- 
diodes. An observation unit consisted of a set of 1000 scans 
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traversing 140 in optical delay digitised into 1024 steps 
of 0.2 ms samples. 

During post-processing of the data, the two signals were 
differenced to mitigate noise introduced by scintiUation and 
the squared visibility was estimated after bias subtraction 
for each observation. A nonlinear functio n was ap plied to 
partially correct for residual seeing effects (|lreland| [2006). 

Reference stars close 7^ Vel on the sky, chosen to pro- 
vide calibration of the system response, were interleaved 
with the target observations. The system response as a func- 
tion of time was quantified using the adopted stellar param- 
eters of the calibrator stars given in Table [T] Measurements 
of were finally produced by linearly interpolating the sys- 
tem response to the time of target observation and scaling 
the observed squared visibility appropriately. This proce- 
dure resulted in a total of 278 estimations of as sum- 
marised in Table [2] The majority of the observations were 
completed in excellent conditions and the atmospheric cor- 
rection to the final values was small with a mean value 
of the order 2 percent. 



3 ORBITAL SOLUTION 

The theoretical response of a two aperture interferometer 
to the combined light of a binary star can be given by 
l|Hanburv Brown et al.lll970l ') 



^2 ^ Vi + + 2/3IV1IIV2I cos(27rb ■ p/A) 



(1) 



where Vi , V2 are the visibilities of the primary and secondary 
respectively and /3 < 1 is the brightness ratio of the two stars 
in the observed bandwidth. The angular separation vector 
of the secondary with respect to the primary is given by p 
(measured East from North), b is the baseline vector pro- 
jected onto the plane of the sky and A is the centre observing 
wavelength. The observed will vary throughout the night 
due to Earth rotation of b and the orbital motion of the bi- 
nary. The Keplerian orbit of a binary star, i.e. p as a function 
of time, can be parameterized with seven elements: the pe- 
riod P, semi-major axis a, eccentricity e, epoch of periastron 
To, the longitude of periastron uj, the longitude of ascending 
node f2 and the inclination i. When using two-aperture opti- 
cal interferometry, the phase of the complex visibility is lost 
and hence, uj and have an ambiguity of 180°. This stems 
from the fact that the identity of the components cannot be 
determined. Radial velocity measurements can be used to 
remove the ambiguity of u but that of Q, remains. Interfer- 
ometers with three or more apertures can use closure phase 
techniques to identify the components and hence remove the 
180°ambiguity of fl. 



3.1 Component Visibilities 

In the simplest case, stars can be modeled by a disc of uni- 
form irradiance with angular diameter 6. The visibility is 
then given by 



V 



2Ji(7r|b|6>/A) 
7\b\e/\ ' 



(2) 



ter are needed to find an estimate of the true angular di- 
ameter. In the case of a compact atmosphere, the correc- 
tions are small but this may not be the case with extreme- 
l imb darkened stars or stars with extended atmospheres 
(|Tango fc David I2OO2I ). In the CEise of emission line stars, 
care must be taken to ensure that only the continuum is ob- 
served (i.e. free of line emission) because a layer that forms 
an emission line cannot be assumed to be at the same radius 
as the continuum. The optically thick stellar wind of WR 
stars not only produces e mission lines but also obscu res the 
hydrostatic core surface (| Moffat fc Marchenkoll 19961 ). 

Over SUSI's wide observing band, emission lines con- 
tributed a significant portion of the detected light. The ex- 
pected spectral response of SU SI (the combination o f the 
700 nm fiher profile and fig. 5 of lDe Marco et al.ll200dl ') was 
used to estimate that four emission lines contribute approxi- 
mately 32 percent of the light received from the WR compo- 
nent during our observations. Three lines are blends of Hel- 
CIII, CIII-Hel and CII-CIII at AA674.1, 706.6, 723.0 nm re- 
spectively. The remaining line is due to Hell at A656.0 nm. 

As each emission line corresponds to a different exci- 
tation/ionisation, it follows that each line could be formed 
at a different radius in the stellar wind . Moreover, emis- 
sion lines may form over a range of radii. iHiUieil (|1989D and 
iDessart et al.l (|2000l ) have used WR models to predict the 
line formation stratifica tion in the stellar wind . Niedzielskil 
( 1994 ) and Schulte- La dbeck. Eenens fc DavQi (|l995l ') have 
shown that observed WR emission lines have differing equiv- 
alent widths and concluded that these differences were evi- 
dence of excitation/ionisation stratifi cation within the WR 
stellar wind. In the case of 'y^ Vel, iHanburv Brown et al.l 
(197Q.) estimated the uniform disc angular diameter of the 
A465.0nm CIII emission region to be 2.05 ± 0.19 mas and 
the analysis of the VLTI included modeling of the relevant 
WR IR emission line layers as uniform discs, e ach with a 
small difference in diameter (|Millour et al.ll2006l ). 

Unfortunately, the radii of the detected emission line 
forming layers have not been measured over SUSI's opti- 
cal range and SUSI's red system presently has no way of 
isolating a single emission line. One possibility is that the 
four detected emission line forming layers are at essentially 
equal radii allowing a simple single uniform disc model of 
the emission. However, a Gaussian irradiance profile may be 
more applicable if the layers are spread in radii. 

We therefore analyse the interferometric data using two 
models of the WR star, each with a continuum core as a disc 
of uniform irradiance with angular diameter 6c- The emis- 
sion line layers are modeled to be either a single uniform disc 
of angular diameter Ov or a Gaussian shell with angular full- 
width half-maximum of Oq contributing Je of the detected 
irradiance. The theoretical visibilities for these models are: 



^^^ 2Ji(7r|b|gc/A) ^ ^JJi{7v\b\ev/X) 



VwR,G ~ (1 — /e 



7rlb|6'c/A 

2Ji(7r|b|gc/A) 
n\b\ec/X 



■K\b\ev/\ 



+ Je exp 



9'h\b[ 



41n2 A2 



(3) 



(4) 



where Ji is a first order Bessel function. Real stars are limb- 
darkened, therefore corrections to the uniform disc diame- 



The O-star primary is modeled as a uniform disc of angular 
diameter 80 as per equation ([2}. Note that the brightness 
ratio /3 in equation ((T} is usually given in reference to the 
continuum irradiance detected in the observing band i.e. free 
from emission lines. As our observing band is contaminated 
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Table 2. Summary of observational data. The night of the observation is given in cols 1 and 2 as a 
calendar date and a mean MJD. Col 3 is the mean orbital phase calculated from the values in Table [3] 
The baseline and the mean projected baseline (in units of metres) are given in Col 4-5 respectively. 
Reference stars and the the number of squared visibility measures for a night are listed in the last two 
columns. 



Date MJD Phase Nominal Projected Reference # 

Baseline Baseline Stars 



2005 Mar 10 


53439.52 


0.26 


80 


75.52 


c 


Pup, 


X 


Car, J Pup, a Pyx 


10 


2005 Mar 11 


53440.48 


0.28 


80 


75.93 


c 


Pup, 


X 


Car, a Pyx 


15 


2005 Mar 12 


53441.48 


0.29 


80 


75.78 


c 


Pup, 


X 


Car, a Pyx 


12 


2005 Dec 17 


53721.70 


0.86 


80 


76.19 


c 


Pup, 


X 


Car, a Pyx 


9 


2005 Dec 18 


53722.69 


0.87 


80 


76.18 


c 


Pup, 


X 


Car, J Pup 


13 


ofinK Flcir' 1Q 
zuuo i-JGC ly 


Oo ( Zo. Dc5 


U.oo 


oU 




/- 




X 


V-jar, J r up 


1 9 


2005 Dec 20 


53724.68 


0.89 


80 


76.19 


c 


Pup, 


X 


Car, J Pup 


12 


2005 Dec 23 


53727.62 


0.93 


80 


76.02 


c 


Pup, 


X 


Car 


6 


2005 Dec 29 


53733.67 


0.01 


80 


76.02 


c 


Pup, 


X 


Car 


14 


2005 Dec 30 


53734.65 


0.02 


80 


76.04 


c 


Pup, 


X 


Car 


10 


2005 Dec 31 


53735.66 


0.03 


80 


75.97 


c 


Pup, 


X 


Car 


16 


2006 Jan 01 


53736.62 


0.05 


80 


76.22 


c 


Pup, 


X 


Car 


9 


2006 Jan 02 


53737.62 


0.06 


80 


76.04 


c 


Pup, 


X 


Car, J Pup 


15 


2006 Jan 07 


53742.64 


0.12 


80 


76.30 


c 


Pup, 


X 


Car 


9 


2006 Jan 08 


53743.62 


0.14 


80 


76.28 


c 


Pup, 


X 


Car 


9 


2006 Jan 09 


53744.65 


0.15 


80 


76.08 


c 


Pup, 


X 


Car, J Pup 


13 


2006 Feb 08 


53774.61 


0.53 


80 


75.39 


c 


Pup, 


X 


Car 


10 


2006 Feb 10 


53776.60 


0.56 


80 


75.35 


c 


Pup, 


X 


Car 


15 


2006 Feb 11 


53777.57 


0.57 


80 


75.11 


c 


Pup, 


X 


Car 


24 


2006 Feb 13 


53779.56 


0.59 


80 


75.13 


c 


Pup, 


X 


Car 


21 


2006 May 11 


53866.37 


0.70 


5 


4.71 


c 


Pup, 


X 


Car 


5 


2006 May 11 


53866.42 


0.70 


30 


27.47 


c 


Pup, 


X 


Car 


4 


2006 Jun 17 


53903.38 


0.17 


5 


4.30 


c 


Pup, 


X 


Car 


7 


2006 Jun 18 


53904.37 


0.18 


15 


12.95 


c 


Pup, 


X 


Car 


6 



with emission lines it is convenient to introduce the rela- 
tive brightness j3' of the WR component (including emis- 
sion lines) to that of the O-star in the observing band. The 
simple relation 



y2 ^ Vi + I3' V4 + 2l3r('4^)\Vi\\V2\ cos(^) 



/3 = (1 - /e)/3' 



(5) 



can then be used to estimate the usual continuum brightness 
ratio. 



3.2 Fitting Procedure and Uncertainty 
Estimation 

Equation ([ij is strictly only valid for observations over very 
narrow bandwidths. For real detection syst ems wide hand- 
width effects can reduce the observed . I Tango fc Davis! 

present corrections which can be applied to wide 
bandwidth interferometric observations of single stars. In 
the case of binary stars, the way in which the wide band- 
width affects the observed squared visibility {V^) is depen- 
dent on the detection system and subsequent calculations 
(fSod cn 1999). For a scanning detection system (such as the 
one used at SUSI) the modulating term in equation ^ will 
be reduced by a factor that is dependent on the separa- 
tion of the component phase centres in delay space and the 
spectral response of the interferometer. Approximating the 
spectral response as a Gaussian of centre wavelength Ao with 
full-width half-maximum AA and for convenience defining 
ip = 2-Kb ■ p/Ao, then the equivalent to equation ^ giving 
V'^ for a binary star for the case of a wide spectral band- 
width is: 



where 



r(i/') = exp 



32 In 2 



(6) 



(7) 



The term r{il>) corresponds to the autocorrelation of the 
Gaussian envelope of the interference pattern. As AA ap- 
proaches zero, i.e. a narrow bandwidth, then equation ([BJ 
reduces to equation 

In the case of 7^ Vel, the emission lines (in an otherwise 
near-Gaussian spectral res ponse determined by t he 700 nm 
filter profile and fig. 5 of iDe Marco et al.l boOOl ) alter the 
shape of the modulated interference fringes, further compli- 
cating the estimation of . Detailed numerical simulations 
showed that the error introduced when approximating the 
spectral response as Gaussian was smaller than the associ- 
ated measurement error. Hence a Gaussian approximation 
of the spectral response was deemed adequate. 

Initial values of i and Q. were found via a coarse search 
of parameter space with the remaining orbital parame- 
ters limited to within 3 st a ndard deviations of the val- 
ues given in ISchmutz et al.l |l993) and the stellar model 
paramete rs {9o, dg, dg / dv, P , Je) were estimated us- 
ing Hanburv Brown et all lligTOl l. lSchmutz et all (|l997l ) and 
lOe Marco et all (|2000l ) 

The final estimation of parameters was completed us- 
ing minimization as implemented by the Levenberg- 
Marquardt method to fit equation ([6]) to the observed V^. 
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As the radial velocity measures have occurred over a long 
time frame of many periods, the spectroscopically deter- 
mined period will be of higher accuracy than that de- 
termined by in t erfero metry. Therefore the period given in 
ISchmutz et alj |l993) was adopted and the remaining or- 
bital parameters were allowed to vary. The angular diam- 
eters were initially fixed to 0.44 mas for the primary (from 
iHanburv Brown et al.|[l970l ) and 0.22 mas for the WR con- 
tinuum core i.e. a continuurrQ radius of 6Rq at the Hip- 
parcos distance of 258 pc. The emission line contribution for 
each model was fixed to 32 percent of the received WR star 
irradiance and the brightness ratio /3' was a free parameter. 
These values are not critical as the resolution of the SUSI 
configuration during observations is such that the WR con- 
tinuum is essentially unresolved. Moreover Vi, V2 & /3' (and 
hence ^o, ^c, do/Ov & -^^e) are coupled and only affect the 
estimated orbital solution from equation (|6]) if highly inap- 
propriate values are adopted. 



When finding the minimum of the manifold, the non- 
linear fitting program calculates the inverse of the covariance 
matrix. The diagonal elements are used to derive the formal 
uncertainties of the fitted parameters. As equation © is 
non-linear and the visibility measurement errors may not 
strictly conform to a normal distribution, the formal uncer- 
tainties may be underestimates. To confirm the accuracy of 
the values derived from the covariance matrix, three uncer- 
tainty estimation methods were adopted. 



Firstly, the model visibilities of each observation were 
subjected to Monte Carlo realisation of the measurement er- 
ror distribution to produce synthetic data sets. Secondly, the 
bootstrap method (jPress et al.lll993 ) was employed to pro- 
duce a second population of synthetic data. By randomly 
sampling with replacement, a synthetic data set is formed. 
The bootstrap method has the advantage that the measure- 
ment error distribution is not assumed to be known as per 
Monte Carlo but may not produce a representative sample. 
The synthetic data sets from each method were used to esti- 
mate the model parameters with the same non-linear fitting 
program, thus building a distribution of each model param- 
eter. 



The final uncertainty estimation method involves a like- 
lihood based random walk through parameter space using 
a Markov chain Monte Carlo (MCMC) simulation imple- 
ment ed with a Metropolis-Hastings algorithm (see Chapter 
12 of i Gregory : 2005a for an introduction to MCMC). This ap- 
proach yields the full marginal posterior probability density 
function (PDF) but may fail to fully explore pathologically 
narro w probability pe aks in a reasonable number of itera- 
tions (lGregorvl l2005bl ) . Furthermore, the current knowledge 
of the system can be included in the analysis by assuming 
an a priori distribution. For example, the effect of adopting 
the spectroscopic period (and associated uncertainty) on the 
remaining orbital parameters can be included into the un- 
certainty estimation of the remaining model parameters. 



Table 3. Comparison of the fitted parameters of this work with 
those found in the literature. 



Parameter 


Unit 


This Work 


Literature 


Ref 


P 


days 


(78.53 ± 0.01)" 


78.53 ± 0.01 


S97 


a" 


mas 


3.57 ±0.05 


4.3 ±0.5 


HB70 


e 




0.334 ± 0.003 


0.326 ±0.01 


S97 


To 


MJD 


50120.4 ±0.4 


50120 ±2 


S97 


LU 


dog 


67.4 ±0.5 


68 ±2 


S97 


n 


deg 
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" adopted parameters, definition is ambiguous, see text for more 

details, '= square brackets deno te a ran ge. 

HBTO: iHanburv Brown et al ] lll970h : SL87: ISt.-Louis et aP 
lll987l) : S97: ISchmutz et al.l lll997l '): DMOOi lDe Marco et al.l (I2OO0I ') 



3.3 Results 

Initial analysis resulted in orbital parameters that, when 
combined with the spectroscopic values as shown in Sec- 
tion |4]2l produced a distance of approximately 330 pc. Con- 
sequently, the angular diameter of the WR continuum core 
was revised to 0.17 mas (consistent with radius expectations 
at this new distance) which had a negligible small effect the 
fitted orbital parameters. We note that the WR core is es- 
sentially unresolved and therefore any diameter uncertainty 
(such as limb-darkening corrections to the equivalent uni- 
form disc angular diameter) are insignificant. 

The results of the two models (uniform disc and Gaus- 
sian shell) were completely consistent in all orbital parame- 
ters and the final values of the fitted parameters are given in 
Table |3] Example data from two nights with the fitted mod- 
els is shown in Fig. [T] and the projected orbit on the plane 
of the sky is shown in Fig. O The reduced of the two fits 
were approximately 2.2 implying that the formal measure- 
ment errors calculated by the data reduction software are 
underestimated by 49 percent. Therefore the measurement 
uncertainties were scaled by 1.49 before the parameter un- 
certainties were estimated using the techniques outlined in 
Section [3.21 The Monte Carlo and bootstrap methods were 
set to each generate 10^ synthetic data sets while the MCMC 
simulations completed 10^ iterations. All combinations of 
model and uncertainty estimation methods produced simi- 
lar distributions of the free parameters which were, for the 
orbital elements, Gaussian in appearance. Furthermore, the 
standard deviation of each model parameter from all un- 
certainty methods were consistent to within rounding. The 
standard deviations derived from the MCMC simulations 
were only slightly larger as the probability space of all pa- 
rameters was explored: the period, WR continuum core di- 

^ The continuum radius corresponding to a mean Ros seland opti- 
cal depth of approximately unity is given in lSchmutz et al. ( 19971 ) 
to be about 6 Rq . 
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Figure 1. Data from the nights of 2005 Mar 11 and 2006 Feb 
13 where each data point represents a measure of with the 
associated formal error. The values of Table |3] have been been 
used to show the uniform disc (solid) and Gaussian shell (dashed) 
emission layer models. 



ameter and emission line contribution varied within a Gaus- 
sian likelihood with assumed values given in Table |3l The 
distributions of the emission region parameters were some- 
what asymmetric with a larger tail at higher values. This is a 
direct result of the small number of short baseline measure- 
ments combined with the coupling of the component visibil- 
ities and the brightness ratios. All uncertainties quoted in 
Table|3]are the standard deviation values of the MCMC sim- 
ulations as we believe they represent the most realistic and 
conservative parameter uncertainty estimates for our data 
set. 

There is no evidence that the visual companion 7^ Ve l or 
the nearby K4 V star detected bv lTokovinin et al] (|l999l ) af- 
fect our data. 7^ Vel is approximately 2.5 magnitudes fainter 
and separated by 41" from 7^ Vel; the K4 V star has a 
separation from 7^ Vel of approxi mately 4.7" and is ap - 
proximately 14.8 magnitudes fainter (|Tokovinin et al.lll999l l. 
These stars are either out of the field-of-view or too faint to 
be detected by SUSI. 

The or bital elements derived from radial velocity mea- 
surements jSchmutz et al.l Il997f ) ar e consistent wi t h our 
data as is the inclina tion inferred bv ISchmutz et al] (Il997l ) 
and iDe Marco et al.l ImOOiV The discrepa nt angular semi- 
major axis of Hanburv Brown et al.l l|l97Ci ) is largely due to 
their adopti on of the inclination (o and i axe co-dependent 
as shown in lHummel et al]|l993l ). This correlation between 
a and i was confirmed by the ad-hoc fixing of the inclination 



Figure 2. The relative orbit (solid line) of the WR-star about 
the O-star projected on the plane of the sky at an orbital phase of 
0.4. The dotted line is the line-of-nodes, the open triangle signifies 
periastron and the grey squares are parts of the orbit that were 
observed with SUSI. The stellar components (filled circles) are 
drawn to scale using the values of Table |3] The full-width half- 
maximum of the Gaussian shell and the angular diameter of the 
uniform disc emission line layer models are shown as the grey 
circles about the WR component. 



as per lHanburv Brown et all ()l970l ) and the resultant angu- 
lar semi-major axis was found to be consistent with their 

value. 

The single angular separation vector of iMillour et al.l 

(|2006D could be use d to remove the 180 ° ambiguity in our 
value of Q. However, iMillour et al.l (I2OO6I ) have stated their 
belief that they have identified the components correctly 
but are cautiously awaiting the final analysis of calibration 
data. Therefore, we tentatively set the position angle of the 
ascending node to that given in Table O Using the final 
orbital values at the epoch of the VLTI measurement, the 
separation of the components is 3.85 ± 0.06 mas and the 
WR component is located 80±1° (measured Eastwards from 
North) relative to the O-star. These values are marginally 
consistent with th e VLTI measur ement of 3.62]'^q'3q mas and 
73tli deg given in IMillour et all (|2006i l . 



3.4 Polarimetry Analysis 

Polarimetric analyses of 7^ Vel have produced the only 
pre-existing estimations of the position angle of the as- 
cending node. However, the parameter, Q, quoted in the 
polarimetry literature is not the position angle of the as- 
cending node but rather the position angle of the rotation 
axis (orbit normal) projected o n to the plane of the sky 
(|Brown. McLean fc Emslidlll978l ). This angle is included in 
the polarimetry analysis to account for the rotation of the 
binary system in relation to the coordinate system of the ob- 
served Stokes parameters Q and U. In the plane of the sky, 
the relationship between the two definitions of f2 is a simple 
difference of 90° . However, the polarimetry ro tation param- 
eter Q is usually defined in the Q-U plane (see [Robert et al.l 
Il989f) and hence an additional factor of two is introduced 
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±90 



(8) 



where Qan and Q,ra are the position angles of the ascending 
node in the plane of the sky and the projected rotation axis 
in the Q-U plane respectively. The addition or subtraction 
of 90° is dependent on the orientation of the binary orbit 
with respect to the North C elest ial pole. Unfortunatel y nei- 
ther 1st .-Louis et all |l983) nor ISchmutz eral] l| 19971 ) give 
details of which p lane their quoted ftra is measured in. For 
the range given in lSt.-Louis et al.l (|l987l ) Qan could be 190°- 
250°or 230°-260°if transformed from the plane of the sky 
or from the Q - U pl ane respectively. Thus, the values of 
ISt.-Louis et~an (|l987l ) are compatible with the interferomet- 
ric value regardle ss of which plane their value is measured. 
The Qan value of lSchmutz et all (|l997l ') could be 232° (plane 
of the sky) or 251° {Q-U plane) respectively. Thus flan de- 
rived from the Q-U plane is close to our value. However, 
without an estimate of the uncertainty we cannot determine 
its consistency with our value. 

As polarimetric data of 7^ Vel are the only other avail- 
able experimental values to estimate the position angle of the 
ascending node and the inclination (a nd to u pdate the po- 
larimetric mass-loss rate determined bvlSt.-Louis et al.lll988l 



in Section l477)) the polarimetry data of lSt.-Louis et al" 
were reanalysed using the model of iBrown et al.l 



1987 ) 



19821 ) 



(with the corrections made bv lSimmons fc Bovlelll984 ) and 
the rotation angle il™ included in the Q-U plane as per 
[Robert et al.l(|l989l ). A minimization was performed with 
the period and eccentricity fixed to the interferometric val- 
ues given in Table [3] The starting values for the remain- 
ing parameters were found using the interferometric values 
given in Tableg] The best fit values were To = 50120.9 MJD, 
Xp = 163°, D.ra = 303°, i = 66°, r, = 0.036, Qo = 0.077 and 
Uo = —0.077. Converting to usual binary star parameters, 
Xp — Lo -\-90° and equation ^ give lo = 73° and Qan = 242°. 
A 10^ iteration MCMC simulation was performed to investi- 
gate the effect that the assumed parameter and eccentricity 
values have on those that were fitted. All parameters with 
the exception of the inclination produced Gaussian PDFs 
centred roughly on the best fit values. However, the PDF 
of the inclination had multiple peaks with 90 percent of the 
probability conta ined between 64 . 5° and 77.5°. This confirms 
the statement bv lSchmutz et al] l|l997h that the inclination 
cannot be accurately constrained by the available polarime- 
try data. The PDF of Qan produces an uncertainty of ap- 
proximately 4° and hence our polarimetric value of Qan is 
consistent with the interferometrically determined value. 



4 DISCUSSION 

4.1 Wind- Wind Interactions 

The orbital values of Table |3] indicate that the angular 
separation at periastron is approximately 1.29 mas sug- 
gesting that the wind of the O-star primary may in- 
teract with the emission line forming layers of the WR 
wind. Moreover, a higher temperature and density col- 
liding winds zone (as inferred by th e X-ray spectrum of 
-Y^ Vel IWiUis. Schild fc Stevenj Il995l : IStevens et al] 1 19961 : 
iHenlev. Stevens fc Pittardll2005l ) may also affect the appear- 
ance of the emission line forming layers. However, the angu- 
lar diameter of these layers, whether modeled as a uniform 



disc or Gaussian irradiance distribution, is not adequately 
constrained by our data. The effect of an irradiance asym- 
metry due to wind-wind interactions on the orbital param- 
eters is considered to be negligibly small as the WR wind is 
expected to dominate that of the O-star. Furthermore, the 
short baseline data, which should be most strongly affected, 
are limited and were collected during an orbital phase when 
the O-star was outside the estimated emission line forming 
layers. 



4.2 Distance 

The distance to a double-lined spectroscopic binary star that 
has a well-determined visual or inter ferometric or bit can be 
found using the dynamical parallax l|Heintzlll978l ): 



(ai -f 02) 



(9) 



The semi-major axis of the primary (secondary) component 
orbit about the system centre-of-mass, ai(a2), is measured 
in astronomical units. 

Using the semi-amplitude of th e component radial ve- 
locity given in lSchmutz et al.l l|l997l ). {Ki = 38.4±2kms"\ 
K2 = 122±2kms~^) and the interferometric orbital values, 
oi and 02 were calculated (in km) using 



ai,2 



43200is:i,2P/r^ 



(10) 



These values are given in Table |4l The resulting dynamical 
parallax is tt^ — 2.97 ± 0.07 mas which yields a distance of 
336t^pc. 

Compared to t he Hi pparcos (jSchaerer et al.l Il997l : 
Ivan der Hucht et all Il997h and the VLTI distance 
(Mill our et al.l 20061 ) our uncertainty in the measured 
distance is lower by an about order of magnitude. More- 
over, there is considerable contrast between the distance 
determinations. It has been shown that the dynamical 
parallax of some binary sta rs can be significantly different 
to the Hipparcos parallax l|Davis et al.n2005l : iTango et all 
[2OOQ) . Binary motion can influence the Hipparcos value and 
where poss i ble, h as been included in the data reduction 
(|Lindegrenl 19971). Th e 7^^ Vel entry in the Hipparcos 
catalogue (|ESA| (l997l ) indicates that it was not treated 
as binary but it is entered into the Variability Annex: 
Unsolved variables. Furthermore, a critical reassessment of 
the Hipparc os catalogue has id entified a number of defects 
in the data (Ivan Leeuwenll2005h and a new reduction of the 
raw data is in progress (|van Leeuwen fc Fan tino 2005,). The 
preliminary revision to the Hipparcos parallax of ■y'^ Vel 
is 3.35 ± 0.29 mas (van Leeuwen, private communication) 
i.e. a distance of 299^24 PC which is marginally consistent 
with our value. The VLTI distance (iMillour et al.l |2006| ) 
relies on a single angular separation measure and assumed 
orbital elements from the literature which are inferior to 
our interferometric values. 

The 7^ Vel Hipparcos parallax raised the question of 
the assumed membership of the Vela 0B2 association due 
to the disparity of their distances: 258 pc vs 450 pc. The Hip- 
parcos proper motion and parallaxes of nea rby OB associa- 
tions including Vela OB2 were analysed bv lde Zeeuw et al] 
(1 19991 ) who concluded that 7^ Vel was an edge member 



with probability of 64 percent. The 93 members produce 
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a mean distance of 410 ± 12 pc (|de Zeeuw et al.lll999l ) i.e. 
a mean parallax of ap proximately 2.44 mas. Using fig. 11 of 
Ide Zeeuw et~al] 1 19991) it was estimated that 60 percent of 
the Vela QB2 members lie within 0.5 mas of the mean value. 
Furthermore. IPozzo et all (|2000| ) model the Hipparcos paral- 
lax dispersion as a Gaussian of standard deviation 0.68 mas. 
Our parallax moves 7^ Vel from the edge of Vela 0B2 to 
well within the association . Using the Galactic coordinates 
017^ Vel in lvan der Huchtl (fioOli ) (and a Sun-Galactic centre 
distance of 8kpc) the revised Galactocentric distances are 
J?GC = 7.97 kpc and zcc = — 45pc. We estimate an uncer- 
tainty in Rgc of 0.01 kpc and zgc of 1 pc which are a factor 
of four less than those estimated with the Hipparcos dis- 
tance. Compared to values found in the lit erature {Rgc = 8 
kpc, ZGC = —35 PC [van der Huchtl I2OOII ) our values place 
7^ Vel farther from Galactic plane but at approximately the 
same distance from the Galactic centre. 



4.3 Component Magnitudes 

The absolute visual magnitude of the system and compo- 
nent stars can be redetermined using our distanc e. There 
is th e possibility of broadband photometry ( Joh nson et al.l 
1 19661 : ICousinsI 1 19721) suffering contamination by the WR 
emission hues (|Smitbl Il968bl : iDe Marco et al.1 mod ) and 
studies of galactic WR distributions and intrinsic param- 
eters almost exclusively use narrowband 'line-free' photom- 
etry (e.g. [Smith 1968bl ; IConti et al.l 1 19831 : IConti fc Vaccal 
|1990D . The mean and standard deviation of narrowband 
apparent magnitude and extinction values found in the 
litera t ure are v = 1.73 ± 0.03 ma g ( van der Hucht et al.l 
1988; Smith, Shara & MoflFat 1990; 'Conti fc Vaccal Il990l : 
van der Hucht_c ^ al. 1997; Do MarcQ et al. 2000l) and A„ = 
0.09 ±0.06 mag Icoiiti fc Vaccalll990l : iDeMarco et al.ll2000l : 
Ivan der Hucht] l200lD . Therefore, combining these values 



with our distance we obtain the system absolute magnitude 
of M„(WR-HO) = -5.99. 

The component absolute magnitudes in both narrow- 
band (w) and broadband (V) can then be found using the 
component brightnes s ratio in v and Af„ — Mv ~ —0.1 mag 
(|van der Huchtl [2001] ). The output brightness ratio /3' from 
the fitting procedure and the estimated emission line con- 
tribution 7e produce a continuum brightness ratio in our 
observing band of /3 = 0.31 ± 0.03. This is consistent with 
the expected brightness ratio at 700nm of 0.3 5 ± 0.04 de- 
termined from fig. 1 of iDe Marco et al. (I2OO0I). We there- 
fore adopt the v brightness ratio of De Marco et al.l (|2000l ) 
rather than scaling the continuum brightness ratio in our ob- 
serving band. We obtain M.„(WR)= -4.33 ± 0.17mag and 
Mv{0)= -5.63±0.10 mag (from A/„(0)= -5.73±0.10 mag) 
where the estimates of the magnitude uncertainties include 
the distance, reddening and the component brightness ratio. 
The uncertainties in previous estimates of the distance to 
7^^ Vel have dominated err or determinations of the c ompo- 
nent absolute magnitudes (Ivan der Hucht et al.lll997l ). This 
limitation has now been removed with our measurement. 



4.4 O-star Component 

The most recent classification of the O - star s pectral type 
is 07.5 given by iDe Marco fc Schmutj 1(199^ . They also 



deduce a bright-giant luminosity class (H) from analy- 
sis of their synthetic Hell line A468.6 nm as other lines 
are unavailable due to assumed metallicity in their model 
and blending effects prevent an observational determina- 
ti on. The recent sprea d-spectrum interferometric analysis 
of iMillour et al.l (|2006l ) has found the O-star models that 
best fit their data have gravity values more in line with 
a super-giant classification. However, the Hipparcos dis- 
tance produces A/v(0)~ —5.1 mag which is more consis- 
tent with a giant or dwarf hence the 07.5 HI-V classifi- 
cation in the VIRh catalogue of galactic Wolf-Rayet stars 
(|van der Hucht|[200ll ). Moreover, the calibration of physical 
parameters of 0-stars has evolved as stellar models develop 
and clear separation into only three luminosity classes I, 
HI & V has become standard. Assuming a spectral type of 
07.5, our Mv(0) is co nsistent with a giant cl assification 
using iHowarth fc Prinial (1989) or lVacca. Garman v & ShulJ 



(|l99d ) but fal ls between a giant and a super-gia nt with the 
calibration of iMartins. Schaerer fc Hillierl (|2005l ). Therefore 
a luminosity class of H-HI for the O-star seems appropriate 
until the knowledge of this area of the HR diagram improves. 

The radius of the O-star is obtained from the measured 
angular diameter and distance to the system was found to 
be 17±2Rq. Th is is in excellent agreement w ith the earlier 
measurement bv iHanburv Brown et al.l (|l970l ). A radius of 
17 Rq falls between all gian t and sup er-giant calibrations 
found in the literature (Howarth fc Prini a 1989: Va cca et al.l 
Il996l : iMartins et al.1 [2OO5I ) supporting a reclassification of 
the O-star to 07.5 H-HI. 

It is more appropriate to use an upper and lower esti- 
mate of the bolo metric cor r ection to find the O-star lumi- 
nosity. Using iMartins et al.l (|2005h to place an upper (lower) 
limit to the bolometric correction (BC) of -3.4(-3.1), we ob- 
tain logL(0)/LQ = 5.51(5.39). Using a mean luminosity of 
logL(0)/I/Q = 5.45 the age of the O-star is estimated at 
(3 .5 ± 0.4) X 10^ yr from the single-star evolutionary models 
of lMevnet et al.1 (|l994h (assum ing Z=0.02) and an e fi'ective 
temperature of T^ff = 35 OOOK (|De Marco et aLlbOOOD . How- 
ever, rejuvenating mass transfer (with WR-star as donor) 
may have occurred in an earlier interactive p hase. Although 
in go od ag reement with the estimates of ISchaerer et al.l 
(| 19971 ) and iDe Marco fc Schmutj (| 19991 ). this age should 
be considered a lower estimate. The age that we have es- 
timated for the O-star and our re-confirmed membership 
of the Vela OB2 association implies that the nearby pop- 
ulation of low-ma ss, pre-main sequence stars detected by 
iPozzo et ai] (|200d ) is coeval with 7^^ Vel. This population h as 
an age range estimated to be 2-6 Myr (|Pozzo et al.ll200ol ). 



4.5 WR-star Component 

The expected M„ of a WC8 st ar has prev iously been 
given as -4.8 to -6.2 mag (|Smithi , 1968b: Co nti et al.lll983l: 
Ivan der Hucht et al.lll988l 'l but more recentlv ivan der Huchtl 



( 20011 ) lias estimated —3.74 mag with a standard deviation of 
0.5 mag. However, only three Galactic WC8 stars have a dis- 
tance of sufficient certainty to be used in the current estima- 
tion of A/„ and previous studies relied almost exclusively on 
j'^ Vel. We therefore individually compare our value to the 
two other Galactic WC8 stars with distance determinations. 
WR135 has = -4.24 and WR113 has M„ = -3.68 but is 
a binary star with persistent dust formation (|van der Huchtl 
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I2001I ). Our value is close to that of WR135 and, as amor- 
phous carbon dust emission is no t seen in the 7^ Vel spec- 
trum (jvan der Hucht et al.l[l99^ . the difference from that 
of WR113 is not significant. 

Bolometric correction f or a W C8 star is given in 
ISmith. Mevnet fc MermiUiodI (Il994l ') as -4.5 derived from 
evolutionary models and cluste r membership. R e cent m odel 
BCs for 7^ Vel are given in iDe Marco et al.l (|200Cl ) and 
range from -3.5 to -4.1. The atmospheric model analysis 
of the galactic WC8 single star WR135 by iDessart et al.l 
l|2000l) produces a BC = -4.0. Due to the non-binary na- 
ture of WR135 and its similarity to the WC8 component 
of 7^^ Vel, we adopt a BC of -4.0. The resultant luminosity 
is logL(WR)/I/Q ~ 5.23. This val ue is in complete accor d 
with the luminosity determined bv iDe Marco et al] (|2000l ). 



4.6 Stellar Masses 

The component masses of a binary star can be extracted 
from the orbital elements using Kepler's third law and the 
ratio of the component semi-major axes about the centre- 
of- gravity: 



(ai -I- 02) 



Ml 
M2 



02 
ai 



(11) 
(12) 



When the semi-major axes are given in astronomical units 
and the period in years the resultant masses are in solar 
units. Using the values determined in Sections 13.31 and 14.21 
the masses of the 7^ Vel components are M{0) — 28.5 ± 
1.1 Mq and M(WR) = 9.0 ± 0.6 M©. There is good agree- 
ment between our values and those derived from recent spec- 
tral analysis (jPe Marco fc Schmutz 1999; Dc Marco et al.i 
I2OOOI ). 

The mass of the O-star is con sistent with a gian t us- 
ing the 'spectroscopic' calibration of IVacc^ et al.1 lIlQQd) but 
is a li ttle low compared to the values of Howarth fc Prinial 
lll989h. A possible exp lanation f or the higher valu e of 
iHowarth fc Prinial (|l989l ') is given in lVacca et al] (119961 ) who 
note that evolutionary models produce masses systemat- 
ically higher than spectroscopic models (the evolutionary 
model analysis in Section [4.41 als o produced a h i gher m ass: 
34 ± 2Mq). With the tables of iMartins et al.1 l|2005l ') the 
O-star mass is consistent with the 'theoretical' giant mass 
but is slightly high for the 'observed' giant mass. Therefore 
our proposed luminosity classification of II-III is also sup- 
ported by our observed mass and O-star calibrations in the 
literature. 

The mass of the WR star determined with the 
mass-lu minosity r elationship given by equation 3 of 
[Schacrc r fc Maeder (1992) is 9.2 Mq (using the estimated 
logI/(WR)/LQ from Section H3|l which is in good agree- 
ment with our mass value. In the catalogue of lvan der Huchtl 
(|200ll ) there are 6 WC stars that have reliable masses and 
range from 9-16 Mq with a mean of 12 ± 3Mq, consistent 
with our value. 

The mass and M„(WR) we have determined can be 
use d to check our cho ice of BC in Section 14.51 Using fig. 
4 of lSmith et al.l (|l994l ). our mass and A/„ produces a point 
that lies above the -4.5 BC line but is lower than the WC7 



and WC9 points. This implies a bolometric correction be- 
tween -3.5 and -4.5, justifying our choice of -4.0 from the 
atmospheric analysis of WR135. 

4.7 Mass-loss 

Previous observational mass-loss estimates of the WR star 
require adjustment using our distances and revised pa- 
rameters before turning to a comparison with the ex- 
pected values from theor y. The mass-loss derived from radio 
continuum at 4.8 GHz of iLeitherer. Chapman fc Koribalskil 
(1992.) with the w i nd ve locity Vac = 1450 kms~^ from 
lEenens fc Williams! (Il994h and corrected for nonthermal 
emission bv IChapman et al.l (|l999l ) is >frad ^ (3.0 ± 
0.5) X 10~^MQyr~^. This is in agreement with the 
average radio-derived mass-loss rate of WC8-9 stars 
of TVlrad = (2 ± 1) X lO '^ Mfi^y r"^ determined by 
ICappa. Goss fc van der Huchtl (12004 ). The polarimetric 
mass-loss of ISt.-Louis et al.l (|l988l ) has been given more at- 
tention due to the improved orbital parameters from the 
interferometric orbit. As descri bed in Section I 3.4[ w e re- 
analysed the polarimetric data of lSt.-Louis et all (|l987l ) but 
restricted all orbital parameters to the values given in Ta- 
ble [3] to obtain r* = 0.035, Qo = 0.077 and Uo = -0.078. 
This corresponds to a polarisation amplitude near periastron 
of Ap ~ 0.0009. Onc e again using Voo = 1450 kms~^ and 
ISt.-Louis et al.l (|l988h the polarimetric mass-loss is found to 
be Mpo\ ^ (8±3) X 10"*^ Mq yr"^ The disparity of the radio 
and polarimetric measurements can be att ributed to density 
variations or clumping in the WR wind ([Moffat fc Robert 
1 19941 ). Clumping affects the radio mass-loss as it assumes 
a smooth wind and is proportional to the wind density 
squared. On the other hand, polarimetric mass-loss is pro- 
port ional to the density and is therefore inse nsitive to clump- 
ing CNugis. Crowther fc Willis .1998 : Hamann fc Koesterkd 
119981 ). The ratio of the above radio and polarimetric mass- 
loss gives an estimate o f a clumping factor m ~ 3.8 ± 1.8. 
iMoffat fc Roberj 11994) predict a value of > 3 and the 
value for the single WC an alysed by Hamann fc Koesterkg 
( 1998 ) is = m = 4. iSchmutz et al.1 (j 19971 ) re-estimated 
Mpoi ^ 7 X 10~^ Mq yr ~"^ and used the Hipparcos A^i-adio 
of ISchaerer et al.l (|l997l ) to find m in the order of 4. The 
consistency with our value is a fortuitous accident of the 
polarimetric mass-loss changing by a small amount and the 
radio mass-loss remaining unchanged at ~ 3 x 10~^ Mq yr~^ 
- the removal of the nonthermal component in the radio fiux 
negated the revision of the distance. 

Theoretical mass-loss rates for a 9Mq (s mooth stellar 
wind) WR are :^ (1.3 ± 0.4) x 10"^ Mq yr"^ (|Doomlll988l : 
iNugis fc LamerdbOOd ). For the dumping- corrected, chemi- 
cal co mposition dependent relationship of lNueis fc Lamer^ 
(j200(]l ) we also find ~ (1.3 ± 0.4) x lO"'^ Mq yr'^ These 
values are consistent within errors of both our mass-loss es- 
timates. Further com parison can be made to the detailed 
modeling of 7^ Vel by jPe Marco et al.1 (j2000h with the use 
of a clumping volume filling factor /y . Our clumping factor 
m is related to /v by y^V = rn~^ (|Hamann fc Koesterkd 
[199^ i.e. /v ~ 7% and gives .M ~ 8 x lO"*^ Mq yr"^ in 
agreement with our polarimetric mass-loss rate. 

The WR wind performance or transfer efficiency 77 = 
Mvaoc/L represents the ratio of radial momentum in the 
stellar wind to the radiation momentum and is < 1 
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Table 4. Physical parameters of 7^ Velorum. 



Parameter 


Unit 


This Work 


Literature 


Reference 


ai 
0,2 
ai 
02 

distance 


lO'^km 
10^ km 

AU 

AU 
mas 

pc 


43.0 ±2.2 
136.6 ±2.3 
0.287 ±0.015 
0.913 ±0.015 
2.97 ±0.07 
336^7 


3.88 ±0.53 
350+5°, 258+31, 368+13 


SSG97,H97 
HB70, SSG97/H97,M06 


R(0) 
L(0) 
M{0) 
age{0) 


map" 

R0 

L© 
M0 
Myr 


—5.63 ± 0.10 
17 ±2 
2.8 X 10^ 
28.5 ± 1.1 
3.5 ±0.4 


—5.10 ± 0.10" 
17 ±3, 12.4 ± 1.7 
(2.1 ±0.3) X 10^ 
30 ±2 
3.59 ±0.16 


DM00 
HB70, DM99 
DM00 
DM00 
DM99 


M„(WR) 
Mv(WR) 
R(WR) 
L(WR) 
M (WR) 

-Mr-ad 

Mpol 


mag 
mag 
R0 

L© 

M© 
Mgyr-i 
Mgyr-i 


-4.33 ±0.17 
-4.23 ±0.17 
(6 ± 3)" 
1.7 X 10^ 
9.0 ±0.6 
3 X 10~5 
8 X lO^'^ 


-3.84 
-3.76 ±0.20" 
6 

(1.7 ±0.4) X 10^ 
9±2 
3 X 10-5 
7 X lO"** 


DM00 
DM00 

S97 
DM00 
DM00 
SSG97 

S97 



" negle cted distance uncertain ty. 

SSG9 7: ISchaerer et al.l lll997l); H97: Ivan der Hucht et al.l lll997^ : DM00: iDe Marco 



19701) : S97: ISchmut2 et al.l lll997fr 



_ _ . . .. . .. . . ^^^^^^^^^Schmutz 

19991) : DM qO:lDe Marco et aP (120001 ): M06: lMillour et al.l (|2009 ): HB70: iHanburv Brown et al 



for winds that are within the sinRle-scatterinK li mit for 
a ra diatively driven fl ow ([Owocki fc GaylevI Il999l ). How- 
ever, ILucv fc Abbott! (|l993l ) have shown n ~ 10 is pos- 
sible for winds that are radiatively driven with multiple- 
scatterins effects. A typ ical v a.lue of ri ~ 10 , give n in 
lOwocki fc GavlevI (|l999l ) while jPe Marco erall (|2000l ) ap- 
proximate r; ~ 7 for 7^ Vel assuming a 10% clumping volume 
filling factor. The wind performance, using our polarimet- 
ric mass-loss and luminos ity with v^o = 1450 kms~^ from 
lEenens fc WiUiamsl l|l994 ). is T] ~ 3.5; consistent with val- 
ues where multiple-scattering eff ects drive the WR wind and 
mass- loss ||Lucv fc Abbottlll993l ). 



centre of the 2-6 Myr range estimated bv lPozzo et"al] (|2000D 
for the nearby population of low-mass, pre-main sequence 
stars, implying coeval formation with 7^ Vel. 

We have shown that the narrow-band absolute visual 
magnitude and mass of the WC8 component is consistent 
with other galactic WR stars of similar type and theoret- 
ical mass-luminosity relationships. The mass-loss rates de- 
termined by radio and polarimetric measurements have been 
revised to include new information in the literature and our 
orbital parameters. These values allowed an estimation of 
the clumping volume filling factor which is used to show 
agreement between recent atmospheric model analysis and 
our clumping insensitive mass-loss rate. 



5 CONCLUSION 

We present the first complete orbital solution for 7^ Vel 
based on interferometric measurements with SUSI. Emission 
line contamination has been simply modeled and found not 
to affect the orbital solution. 

In combination with the latest radial velocity measure- 
ments, the dynamical parallax and distance to 7^ Vel have 
been revised. The constraints we have placed on the Galac- 
tic location of 7^ Vel are the tightest to date by an order 
of magnitude and membership of 7^^ Vel in the Vela OB2 
association has been confirmed. The contrast between our 
distance and previous measurements in the literature is such 
that all distance-dependent fundamental parameters require 
revision e.g. absolute magnitudes, mass-loss rates. Moreover, 
the formal uncertainty in the distance can now be included 
in the determination of fundamental parameters. 

Using calibrations found in the literature, the radius, 
absolute visual magnitude and mass of the component O- 
star are consistent with a luminosity classification of II-III. 
The age that we have estimated for the O-star falls at the 
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